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摘 要 : 干旱 是 农作物 生长 发 育 的 主要 环境 胁迫 因子 ,也 是 制约 农业 丰产 丰收 的 关键 自然 要 素 。 农 
业 干 旱 监 测 通 常 基于 气象 站 点 观测 数据 ,这 在 一 定 程度 上 难以 反映 区 域 尺 度 的 农业 干旱 状况 。 以 
甘肃 省 为 研究 区 ,基于 MODIS TRMM ESA CCI 等 遥感 数据 产品 和 气象 站 点 数据 ,利用 随机 森林 回 
归 模 型 构建 综合 气象 干旱 指数 (CMDI) ,并 对 甘肃 省 2011 一 2019 年 农作物 生长 季 (4 一 9 月 ) 旱 情 时 
空格 局 及 变化 规律 进行 分 析 。 结 果 表 明 :(1) CMDI 与 实测 值 的 决定 系数 (及 ) 在 各 月 均 高 于 0.634， 
且 与 标准 化 降水 蒸 散发 指数 (SPEI) 在 空间 上 具有 一 定 的 相关 性 ,表明 该 指数 可 反映 农业 干旱 的 发 


生发 展 过 程 ,(2) 甘肃 省 农业 干旱 呈现 明显 的 地 域 分 异 规律 ,干旱 程度 由 东南 向 西北 逐渐 加 重 ,其 


中 河西 地 区 多 为 特 旱 区 和 重 旱 区 , 陇 中 地 区 为 重 (中 ) 旱 区 ,陇南 、 陇 东 、 甘 南 地 区 为 于 旱 - 无 旱 波动 
变化 区 。(3) 2011 一 2019 年 甘肃 省 农业 干旱 在 年 \ 月 太 度 上 均 呈 现 较 大 的 波动 趋势 ,其 中 2012 年 干 
旱 程 度 最 轻 ,2017 年 则 最 为 严重 ;甘肃 省 大 部 分 地 区 在 4 月 和 6 月 , 陇 东 、 陇 南 地 区 分 别 在 5 月 和 9 
月 以 及 甘南 地 区 4 一 9 月 农业 旱情 有 所 减轻 外 ,其 余地 区 在 农作物 生长 季 的 旱情 呈 加 重 趋势 。 
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在 全 球 气候 变 暖 背景 下 ,极端 天 气 事件 频 发 ,其 
中 干旱 因 持 续 时 间 长 .影响 范围 广 而 备 受 关注 一。 
农业 干旱 是 指 在 农作物 生长 过 程 中 , 因 多 种 因素 导 


2S BIOS C CSPED 、 帕 尔 默 干旱 指数 (PDSI) .土壤 水 
分 亏 缺 指数 (SWD1) 、 归 一 化 植被 指数 (NDVI) 降水 
状态 指数 (PCI) .植被 健康 指数 (VHI) 、 植 被 干旱 响 


致 的 水 分 不 足 使 得 农作物 的 正常 需 水 量 受 限 ,从 而 
引发 的 农作物 大 量 减 产 。 干旱 不 仅 对 作物 生产 力 
造成 严重 直接 影响 ,而 且 对 人 类 生计 和 居民 收入 产 
生 间 接 影 响 ,一 直 是 全 球 备 受 关注 的 环境 问题 之 
一 。IPCC AR6 报告 指出 即使 全 球 升 温 稳定 在 1.5~ 
2.0 % ,大 部 分 地 区 仍 将 经 历 严重 且 频 繁 的 农业 生态 
干旱 ”干旱 虽然 是 一 种 无 法 消除 的 天 气 现象 ,但 
若 能 获得 与 作物 状况 有 关 的 实际 时 空 信息 , 则 可 对 
其 进行 有 效 监测 和 管理 从 而 减轻 干旱 影响 "。 
农业 干旱 是 一 种 复杂 的 过 程 ,其 发 生发 展 受 多 
种 不 确定 因素 的 影响 呆 。 随 着 遥感 技术 的 快速 发 


应 指数 (VegDRI) 、 规 模 干旱 条 件 指数 (SDCI) Jh 
温度 (LST) WERE (PET) .土壤 水 分 含量 (SM) 
和 微波 综合 干旱 指数 (MIDI) 等 中 ,一 些 学 者 通过 集 
成 这 些 指 数 以 指示 农业 干旱 状态 ,如 胡 脑 飞 等 " 基 
于 MODIS 和 TRMM 数据 分 析 了 黄土 高 原 的 农业 干 
旱 变化 特征 和 干旱 频率 ,表明 SDCI 指 数 可 有 效 监 测 
该 地 区 作物 生长 季 的 干旱 状况 ; 温 庆 志 等 "构建 了 
综合 遥感 干旱 监测 模型 (IRSDI) ,并 将 其 应 用 于 淮 
河流 域 农业 干旱 时 空 特征 分 析 ; 刘 高 鸣 等 3 利用 决 
策 树 模 型 融合 了 大 气 异常 .植被 状态 .土壤 水 分 等 
言 息 , 相 对 于 单一 指标 在 一 定 程度 上 提高 了 监测 精 


展 ,农业 干旱 监测 逐渐 向 定量 化 ,空间 化 、 多 因子 毕 
合集 成 方向 发 展 *。 目 前 较为 成 熟 的 监测 农业 干 
旱 的 指数 主要 有 标准 化 降水 指数 (SPI) ,标准 化 降水 
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等 的 异常 信息 ;Wu 等 在 VegDRI 指 数 基 础 上 , 提 
出 综合 地 表 干 旱 指 数 (ISDI)。 这 些 农业 干旱 监测 模 
型 与 指数 的 建立 为 开展 不 同 区 域 农业 旱情 监测 与 
评估 提供 了 有 益 借鉴 和 参考 。 

随机 森林 (RF) 是 一 种 基于 分 类 回归 树 的 组 合 
模型 , 较 其 他 算法 具有 处 理 数据 维度 高 .训练 速度 
快 . 抗 干扰 能 力 强 .可 平衡 数据 误差 等 优势 “7”"。 近 
年 来 ,一 些 学 者 将 随机 森林 模型 应 用 于 干旱 预测 研 
究 , 如 沈 润 平等 "9 利用 随机 森林 思想 构建 了 干旱 遥 
感 监测 模型 ; 董 婷 等 ”基于 多 源 遥 感 数 据 , 利 用 随 
机 森林 模型 构建 多 因子 集成 的 旱情 状态 指数 (ID- 
CD ,对 我 国 北方 旱情 状况 评估 表明 该 指数 在 大 面积 
旱情 监测 中 应 用 性 较 强 。 甘 肃 省 大 部 分 地 区 位 于 
干旱 半 干 旱 区 ,生态 环境 脆弱 ,是 农业 干旱 的 高 频 
发 生地 中 。 以 往 对 于 甘肃 省 农业 干旱 的 研究 多 基 
于 气象 站 点 数据 或 单 因子 遥感 指数 ,在 区 域 农 业 旱 
情 监测 与 评 佑 方面 尚 存在 很 大 不 足 。 本 文 应 用 MO- 
DIS TRMM ESA CCI 等 多 源 遥 感 数据 产品 ,结合 气 
伺 ( 降 水 和 LST) 水 文 (土壤 湿度 ) 植被 (植被 状态 ) 
和 地 形 (高 程 坡度 、 坡 向 ) 等 因素 ,在 计算 不 同 干旱 
表征 因子 和 综合 气象 干旱 指数 (CMDI) 基 础 上 ,构建 
甘肃 省 农业 干旱 监测 随机 森林 模型 ,并 对 甘肃 省 
2011 一 2019 年 农业 旱情 时 空格 局 及 变化 规律 进行 
分 析 , 从 而 为 认识 甘肃 省 农业 干旱 发 生发 展 规律 及 
西北 和 干旱 区 农业 和 干旱 预测 提供 方法 支撑 和 科学 
参考 。 


1 研究 区 概况 


HIA (92213' -108?46'E,32?31' -42557' N) Jl 
处 西北 内 陆 , 为 青藏 高 原 .黄土 高 原 和 内 蒙古 高 原 
的 交汇 地 带 "。 甘 肃 省 地 形 呈 狭长 状 ,东西 长 约 
1600 km, 南 北 宽 约 500 km, £ H {E 42.58X 10* km2; 
地 势 自 西南 向 东北 倾斜 ,平均 海拔 1500 m 以 上 ,以 
乌 鞘 岭 为 界 分 为 外 流 区 和 内 流 区 ;气候 类 型 多 样 ， 
从 南 向 北 依次 划分 为 北 亚热带 湿润 区 、 暧 温带 湿润 
区 、 高 山高 寒 区 、 冷 温带 半 湿 润 区 和 冷 温带 半 干 早 
区 5 个 气候 带 下 。 按 照 自然 地 理 条 件 差异 ,甘肃 省 
通常 划分 为 河西 Bn EA DEAA H RS 5 T LOCI 
(图 1)。 在 中 国 九 大 农业 区 中 ,甘肃 省 农业 区 大 部 
分 位 于 北方 干旱 半 和 干旱 区 , 少 部 分 位 于 青藏 高 原 区 
和 黄土 高 原 区 。 根 据 甘 肃 省 第 三 次 全 国 国土 调查 
数据 ,2019 年 甘肃 省 耕地 面积 约 0.32x10' km’, EL 


注 : 该 图 基于 国家 测绘 地 理 信 息 局 标准 地 图 服务 网 站 下 载 的 审 图 
号 为 GS(2020)4630 号 的 标准 地 图 制作 , 底 图 边界 无 修改 。 下 同 。 
图 1 甘肃 省 农业 分 区 


Fig. 1 Agricultural division of Gansu Province 


多 位 于 陇 中 地 区 和 陇 东 地 区 ,以 及 河西 地 区 山 前 绿 
洲 和 陇南 地 区 河谷 地 带 。 


2 数据 与 方法 


2.1 数据 来 源 

采用 1971 一 2019 年 甘肃 省 26 个 气象 站 点 的 月 
降水 量 和 月 均 温 数据 ,其 来 源 于 中 国 气象 数据 网 
(http://data.cma.cn/) ; MODIS 数据 采用 2004—2019 
年 植被 指数 产品 NDVI(MOD13A2,16 d, 1 km) .LST 
产品 (MOD11A2,8 d, 1 km) 和 2015 年 土地 覆盖 类 型 
产品 (MCD12Q1,500m) ,其 中 MCD12Q1 用 于 甘肃 
省 农田 区 的 提取 ,3 种 数据 产品 均 来 源 于 美国 国家 
航空 航天 局 网 站 (https://ladsweb.modaps.eosdis.nasa. 
gov/) ,涉及 行列 号 h25v04、h25v05 FII h26v05 ; ERKE 
水 数据 采用 2004 一 2019 年 TRMM 3B43 月 尺度 数 
据 , 空 间 分 辨 率 为 0.25°x0.25° ,来 源 于 美国 国家 航 
空 航天 局 网 站 (https://disc.sci.gsfc.nasa.gov/); 土壤 湿 
度数 据 采 用 2004 一 2019 年 欧 空 局 “气候 变化 倡议 ” 
Ii H (ESA CCD) 提 供 的 全 球 土壤 湿度 主 、 被 动 微波 
集成 数据 集 , 空 间 分 辩 率 为 0.25°x0.25° ,时 间 分 辨 
率 为 1 d, 该 数据 集 来 源 于 欧 空 局 网 站 (https:/www. 
esa-soilmoisture-cci.org/ ) ;地 形 数据 产品 采用 SRTM 
DEM 数据 ,空间 分 辨 率 为 90 m, 该 数据 来 源 于 地 理 
空间 数据 云 (http://www.gscloud.cn/)。 
2.2 TRMM 与 ESA CCI 数 据 降 尺 度 

考虑 到 甘肃 省 农作物 一 般 为 夏收 和 秋收 ,因此 


202302.00276v1 


chinaXiv 


将 研究 时 段 集 中 于 农作物 生长 季 (4 一 9 月 )。 由 于 
不 同 数据 空间 分 辨 率 差 异 较 大 ,首先 对 研究 区 相应 
时 段 的 MOD13A2、MOD11A2 采 用 平均 值 法 合成 为 
月 尺度 上 的 数据 ,采用 双 线 性 插值 法 将 SRTM DEM 
重 采样 为 1 kmx1 km; 然 后 将 ESA CCI 每 日 湿度 数 
据 合 成 为 月 尺度 数据 ,并 将 空间 分 辨 紊 同 为 0.25°x 
0.25° 的 TRMM 3B43 和 ESA CCI 月 数据 降 尺 度 至 
1 kmX1 km. 

TRMM 数 据 提供 的 三 维 降水 分 布 信息 已 广泛 应 
用 于 国内 降水 研究 中 7, 但 其 空间 分 状 率 较 低 ,小 
区 域 尺度 内 像 元 属性 值 相 同 。 为 了 提高 TRMM 数 
据 空间 分 辩 率 ,根据 李 净 等 所 提出 的 TRMM 数据 降 
尺度 方法 ,以 经 度 .纬度 海拔 .NDVI 等 与 降水 数据 
相关 性 较 强 的 变量 作为 线性 回归 模型 的 自 变量 ,将 
TRMM 3B43 数据 降 尺 度 至 1 km; 同 理 采用 Yu 等 ” 
提出 的 土壤 湿度 数据 方法 ,利用 高 分 辨 率 的 NDVI 
和 LST 将 低 分 辩 率 的 土壤 水 分 进行 分 解 ,并 将 ESA 
CCI 数 据 降 尺 度 至 1 km。 具 体 过 程 如 下 (以 TRMM 
数据 为 例 ): 

P-a,*a,Lon* a,Lat * a,DEM * a,NDVI*& (1) 
式 中 :P 为 降 尺 度 后 的 降水 值 (mm);Lon 和 Lat 分 别 
为 经 度 (°) 和 纬度 (°);DEM 为 海拔 (m);NDVI 为 归 
一 化 植被 指数 ,由 MOD13A2 产 品 中 获取 ;ao~as 为 线 
性 回归 拟 合 系数 ,其 中 a 为 常数 项 ;ge 为 残 差 。 首 先 
T£ DEM 和 NDVI 重 采样 至 低空 间 分 辨 率 0.25°, 预 估 
降水 量 (Pii) 多 元 线性 回归 方程 表示 如 下 : 

P=astaLontaLat+a,DEM+a,NDVI (2) 

KEJK ZR 2E (ATir, mm) W KEJK Œ (Toras , mm) 5 P 
估 降 水 量 (Pa,mm) 之 差 , 计 算 公 式 如 下 : 

AT ir ST Pa (3) 

将 式 (2) 得 到 的 拟 合 系数 代入 至 空间 分 辨 率 为 
1 kmx1 km 的 NDVI 和 DEM 中 ,得 到 高 分 辩 率 的 预 
WREIK EE (Par, mm) ,计算 公式 如 下 : 

P,,7a,*a,Lon*a,Lat- a,DEM * a,NDVI (4) 


将 AT 插值 获得 空间 分 辩 率 为 1 kmx1 km 的 
ATm ,最 后 将 其 与 Pm 相 加 得 到 降 尺 度 后 的 降水 数据 
P: 

P=P n +AT (5) 


23 不 同 干旱 指数 计算 

农业 干旱 不 仅 涉及 土壤 水 分 ,降水 量 作物 的 
生长 状态 及 冠 层 温度 ,还 与 地 貌 类 型 和 土地 覆盖 类 
型 等 因素 有 关 m。 而 植被 状态 指数 (VCI) 是 表征 农 
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作物 生长 状态 的 关键 因子 , 它 对 西北 大 部 分 历史 干 
旱 的 空间 监测 效果 好 ,干旱 空间 演变 特征 显著 ; 温 
度 状态 指数 (TCI) 是 表征 农作物 冠 层 温度 的 关键 因 
子 ,在 干旱 状态 下 , 冠 层 温度 升 高 ,作物 所 需 的 水 分 
减少 , 则 会 影响 作物 的 正常 生长 ,从 而 受到 干旱 胁 
迫 ; 降 水 状态 指数 (PCI) 由 TRMM 3B43 月 降水 数据 
归 一 化 处 理 得 到 ,实现 对 降水 短缺 数据 的 监测 ,可 
以 表征 连续 的 降水 量 异 常 信息 ;土壤 湿度 决定 农 作 
物 的 水 分 供应 情况 ,土壤 干旱 直接 由 土壤 湿度 过 低 
造成 ,由 于 微波 遥感 的 一 系列 优势 ,所 以 将 月 土壤 
湿度 数据 标准 化 处 理 得 到 土壤 湿度 状态 指数 
(CSMCI) 作 为 表征 土壤 湿度 异常 信息 的 关键 因子 ; 综 
合 多 种 遥感 干旱 指数 的 农业 干旱 监测 , 比 单纯 遥感 
因子 具有 更 大 优势 ,温度 植被 干旱 指数 (TVDI) 是 由 
LST-NDVI 构 成 的 特征 空间 且 与 土壤 水 分 具有 强 相 
关 性 ,不 仅 能 够 克服 影像 分 辨 率 小 的 特点 ,可 以 较 
好 地 反映 土壤 湿度 信息 ;微波 综合 干旱 指数 (MIDI) 
综合 PCI、TCI、SMCI 3 个 因子 ,同时 考虑 降水 温度 、 
土壤 湿度 ,是 表征 干旱 的 关键 因子 ;高 程 (DEM) \ 坡 
HE (SLOPE) . 坡 向 (ASPECT) 对 区 域 干旱 同样 有 重要 
影响 。 采 用 张 静 等 "研究 西北 干旱 区 确定 的 权重 
值 来 计算 MIDI( 表 1) ,确定 PCI 的 权重 值 为 0.6,TCI 
的 权重 值 为 0.2,SMCI 的 权重 值 为 0.2。 

CMDI 不 仅仅 考虑 了 降水 和 莹 发 因子 ,而且 考 
虑 了 近 30 4 和 近 90d 的 SPL, 时 间 尺 度 较 好 , 比 利 用 
单 因 子 的 气象 干旱 指数 具有 更 大 优势 ,因此 利用 
CMDI 作 为 模型 的 输出 变量 H。SPEI 作 为 分 析 干 旱 
趋势 的 理想 指标 ,能 表征 一 个 地 区 干 湿 状 态 偏离 常 
年 正常 水 平 的 程度 。 为 避免 研究 时 段 内 突 发 极端 
天 气 事件 对 指数 造成 影响 ,本 研究 根据 月 降水 量 和 
月 均 温 计算 1971 一 2019 年 月 尺度 CMDI 5 SPEI 以 
使 得 各 值 在 正常 范围 内 。 以 CMDI 和 SPEI 的 干旱 等 
级 划分 表征 干旱 程度 ,其 划分 标准 参照 国家 标准 ' 当 
( 表 2), 并 利用 站 点 实测 CMDI 与 SPEI 来 评价 模型 
模拟 结果 。CMDI 指 数 的 计算 方法 如 下 : 

CMDI= aSPL,  bSPL, + cM, (6) 

式 中 :CMDI 为 30d 的 综合 气象 干旱 指数 ;SPL。 SPIw 
为 近 30 d 和 90 d KY SPL; Ma X 30 d 的 相对 湿润 指 
数 。 本 文 根 据 张 调 风 等 研究 成 果 , a 平均 值 取 0.4; b 
平均 值 取 0.4; c 平均 值 取 0.8 富 。 本 文采 用 Vicente- 
Serrano ^ 等 在 SPI 基 础 上 提出 的 SPEI 的 计算 方法 ， 
具体 计算 流程 不 在 此 列举 。 
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表 1 干旱 指数 计算 方法 


Tab.1 Calculation method of drought index 


干旱 指数 计算 公式 


IAN 


植被 状态 指数 (VCI) 


max ”min 


温度 状态 指数 (TCI) 


min) VAY ma 


降水 状态 指数 (PCI) 


土壤 湿度 状态 指数 SMCI=(SM - SM,,)/(SM,.. - SM,,) 

(SMCI) 

温度 植被 干旱 指数 TVDI= (LST - LST") /ELST -IST ) 
(TVDI) 


微波 综合 干旱 指数 MIDI = aPCI + GTCI (1 - a- g)SMCI 


PCI Z(TRMM - TRMM, /(TRMM,,. - TRMM,,,) 


NDVI 为 归 一 化 植被 指数 ;NDVI,,、NDVI, 分 别 为 某 月 归 一 化 植 
被 指数 的 最 大 值 和 最 小 值 。 

LST 为 地 表 温 度 即 地 面 的 温度 (%C);LST,,、LST, 分 别 为 某 月 地 
表 温 度 的 最 高 值 和 最 低 值 (% )。 

TRMM 为 卫星 降水 数据 ,表征 降水 量 的 大 小 (mm:h') ;TRMM、 


SM 为 土壤 湿度 ;SM、SM, 分 别 为 某 月 土壤 湿度 的 最 大 值 和 最 
小 值 。 

LST®Y | pspV" 分 别 为 某 月 NDVI 对 应 的 最 高 和 最 低地 表 温 度 
(C), 

Qa.B 分 别 为 PCI、TCI 的 权重 值 。 


(MIDI) 


表 2 综合 气象 干旱 指数 (CMDI) 及 标准 化 降水 蒸 散 发 


指数 (SPET) 干 旱 等 级 划分 
Tab.2 Classification of drought with respect to CMDI 
and SPEI 
等 级 ”类 型 ” 综合 气象 干旱 指数 ”标准 化 降水 蒸 散 发 指数 
1 无 早 -0.6<CMDI -0.5«SPEI 
2 轻 量 。 -L2«CMDI«-0.6 -].0«SPEI«-0.5 
3 中 旱 -1.8<CMDI<-1.2 -1.5<SPEI<-1.0 
4 重量 -24«CMDI«-1.8 -2.0<SPEI<-1.5 
5 特 旱 ^ CMDI<-2.4 SPEI<-2.0 


2.4 随机 森林 模型 构建 

对 于 非 线性 数据 运用 随机 森林 模型 可 以 减少 
均 方 根 误差 ,取得 较 好 的 拟 合 结果 ”“””。 本 文 以 多 
源 遥 感 数据 提取 干旱 表征 因子 (VCL、TCIL、PCI、 
SMCI, TVDI MIDI, DEM , SLOPE , ASPECT) 为 模型 
输入 变量 ,CMDI 为 模型 输出 变量 ; L} 2004—2010 4E. 
4—9 月 站 点 数据 为 训练 集 ,2011 一 2019 年 4 一 9 月 
甘肃 省 农田 区 为 测试 数据 集 构建 模型 ;利用 2011 一 
2019 年 4 一 9 月 站 点 实测 CMDI 值 与 模型 预测 结 
进行 验证 及 能 力 评价 ,并 以 SPEI 验 证 其 在 空间 上 的 
可 靠 性 ,模型 构建 流程 如 图 2 所 示 。 利 用 python 3.8 
中 Pandas 库 以 及 数据 分 析 工 具 包 Sklearn ,通过 代码 
读 取 训练 集 .测试 集 文本 文件 实现 随机 森林 回归 模 
型 来 估算 CMDI。 


3 结果 与 分 析 


3.1 模型 模拟 结果 评价 
图 3 展示 了 测试 集 输 入 模型 后 站 点 模拟 结果 与 


站 点 实测 CMDI 的 相关 性 。 显 然 ,模拟 结果 与 站 点 
实测 数据 计算 得 到 的 CMDI 在 6 月 的 相关 性 最 高 , 模 
型 拟 合 优 度 达 到 0.935, 其 均 方 根 误差 (RMSE) 和 平 
均 绝 对 误差 (MAE) 均 较 小 ;其 他 月 份 模型 拟 合 优 度 
均 在 0.634 以 上 ,表明 模型 具有 较 强 的 预测 能 

由 于 采用 了 随机 采样 ,训练 出 的 模型 方差 较 小 ,所 
有 月 份 的 RMSE 均 未 超过 0.628 ,说 明 模型 的 泛 化 能 
力 较 好 ,未 出 现 过 拟 合 情况 。 为 验证 模型 预测 值 在 
空间 上 的 可 靠 性 ,将 模型 预测 值 与 同期 的 SPEI 做 相 
关 分 析 ( 图 4) ,在 空间 上 大 部 分 地 区 与 SPEI 具有 很 
好 的 相关 性 ,6 月 的 相关 性 最 高 ,大 致 都 在 0.6 以 上 ; 
其 他 月 份 的 大 部 分 地 区 相关 性 在 0.5 以 上 ,其 中 甘 
南 地 区 每 月 的 相关 性 都 在 0.5 DA. E TRE P rp npe 
南部 、 陇 东部 分 地 区 相关 性 在 0.8 以 上 。 部 分 地 区 
的 相关 系数 为 负 值 且 较 低 ,4.6.9 月 的 负 值 主 要 集 
中 在 河西 中 西部 地 区 ,主要 原因 可 能 是 该 地 区 的 农 
作物 稀少 ,存在 着 很 多 的 不 确定 性 ;5.7 月 的 负 值 主 
要 集中 在 陇 中 东南 部 地 区 ,原因 可 能 是 该 地 区 5 月 
易 发 生 霜 冻 ,7 月 易 发 生 冰 起 ,因此 对 农业 造成 一 定 
危害 ;8 月 的 负 值 集中 在 陇南 东南 部 地 区 ,原因 可 能 
是 该 地 区 人 类 活动 频繁 ,影响 较 大 所 导致 。 结 果 表 
明 ,模型 预测 值 监测 农业 干旱 是 可 靠 的 ,对 于 甘肃 
省 农田 区 域 干旱 具有 一 定 的 监测 能 力 。 

3.2 甘肃 省 农业 干旱 时 空 特征 

3.2.1 农业 干旱 室 间 分 布 格局 及 变化 趋势 ”不 同年 
份 的 甘肃 省 农作物 生长 季 干 旱 空间 分 布 (图 5) 显 
示 , 各 区域 干旱 程度 的 轻重 程度 差异 很 大 。 总 体 
上 ,干旱 高 发 地 主要 集中 在 河西 地 区 和 陇 中 地 区 , 
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实测 气象 数据 


气象 数据 的 获取 


Y 


气象 站 点 月 均 
温 、 月 降水 量 


干旱 指数 的 计算 


F 上 壤 湿 度 状 态 指数 ;TVDI 为 温度 植被 


旱 指数 ;SPEI 为 标准 化 降水 薰 散 发 指数 。 下 同 。 
图 2 随机 森林 模型 构建 流程 图 


Fig.2 Flow chart of random forest model construction 
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注 :MOD13A2 为 植被 指数 产品 归 一 化 植被 指数 (NDVI) ;MOD11A2 为 地 表 温度 (LST) 产 品 ;MCD12Q1 为 2015 年 土地 覆盖 类 型 产品 ;TRMM 
形 数据 ;ESA CCI 为 土壤 湿度 数据 ;DEM .SLOPE .ASPECT 为 高 程 .坡度 . 坡 向 ;VCI 为 植被 状态 指数 ;TCI 为 
F 旱 指数 ;MIDI 为 微波 综合 干旱 指数 ;CMDI 为 综合 气象 干 


3 - (054 3 - (06H 
2 L x-0.949x-0.079 2 上 y-0.972x*0.003 . 
R=0.749 —.. 5 1 R=0935 
1 "mE. S 
0 
0 ô 
rri m zc 
E RMsE-046] € E EE 
s LLL MAE20917 , i 
-3-2-10 1 2 3 -4-3-2-10 12 3 4 
实测 CMDI 实测 CMDI 
4 - (08H 2 - (9H 
- y70.874x-0.042 f y=0.612x—0.081 , 
i : A 2 一 " . 
A R=0.731 , aif3% B R?-0.634 Bu 
0 O 0 
ES ESI "E 
E RMSE-0.5880 SK | S ' RMSE-0.628 
-3 MAE-0.460 MAE-0.493 
4 L L L L L L L f 
-4-3-2-10 123 4 = 2-10 1 2 3 
实测 CMDI 实测 CMDI 


注 :RMSE 为 均 方 根 误差 ;MAE 为 平均 绝对 误差 。 


图 3 4 一 9 月 模拟 CMDI 与 气象 站 点 实测 CMDI 散 点 图 


Fig.3 Scatter diagrams of simulated CMDI and measured CMDI of meteorological stations from April to September 
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图 例 一 省 界 — 农业 区 界 


相关 系数 = -0.1~0.5 


0.5~0.6 = 0.6~0.7 = 0.7~0.8 = 0.8~1.0 


图 4 4 一 9 月 模拟 CMDI 与 SPEI 的 空间 相关 性 
Fig.4 Spatial correlation between simulated CMDI and SPEI from April to September 


甘南 .陇南 、 陇 东 地 区 的 和 干旱 程度 较 轻 。2011 一 
2019 年 ,甘肃 省 农业 和 干旱 程度 最 为 严重 的 年 份 是 
2017 年 ,最 轻 的 年 份 是 2012 年 ;甘肃 省 农业 干旱 呈 
现 较 大 的 波动 趋势 ,年 际 变 化 较为 复杂 ,但 呈现 出 
明显 的 地 域 分 异 规律 , 即 由 甘肃 省 东南 部 向 西北 部 
干旱 程度 加 重 。 这 与 庞 素 非 等 ”基于 TVDI 得 到 的 
甘肃 省 农业 旱 区 旱情 特征 结果 大 体 一 致 。 

从 各 年 份 来 看 ,2011 年 河西 中 部 (张掖 - 嘉 峪 
关 ) 部 分 地 区 和 陇 中 及 陇 东 地 区 是 干旱 最 为 严重 的 
区 域 ,表现 为 重 旱 和 特 旱 ,河西 东部 也 伴随 有 少量 
的 中 旱 现 象 ,陇南 地 区 的 和 干旱 程度 最 轻 ,表现 为 轻 
旱 和 无 旱 。2012 年 除 河 西 地 区 表现 为 重 旱 和 中 旱 
外 ,几乎 甘肃 省 整个 农田 区 都 无 干旱 发 生 , 可见 河 
西 地 区 仍然 是 干旱 发 生 的 高 频 区 域 ; 相 较 于 2011 年 
除 河西 小 部 分 地 区 和 陇南 地 区 干旱 加 重 外 ,其 余地 
区 的 农业 干旱 得 到 明显 的 减缓 , 陇 中 和 陇 东 地 区 没 
有 干旱 发 生 ,陇南 整个 地 区 呈现 干旱 加 重 , 但 表现 
为 轻 旱 。2013 年 河西 东部 干旱 加 重 最 为 强烈 ,表现 
为 特 旱 ,其 余 农田 区 域 均 呈 现 减 轻 趋 势 。2014 年 部 
分 地 区 和 干旱 程度 加 重 , 主要 集中 在 陇南 地 区 ,表现 
为 重 旱 和 中 旱 。2015 年 陇 中 地 区 发 生 了 较为 严重 


的 农业 干旱 ,河西 地 区 的 干旱 有 所 减轻 。2016 年 干 
旱 最 为 严重 的 区 域 是 河西 地 区 的 东部 PRR, fi D 
中 、 陇 南 . 陇 东 地 区 均 有 轻 度 干旱 发 生 。2017 年 整 
个 甘肃 省 农田 区 均 有 干旱 发 生 , 特 旱 、 重 旱 仍然 集 
中 在 河西 地 区 , 轻 旱 和 中 旱 发 生 在 陇南 和 陇 东 地 
区 ;根据 "2018 年 甘肃 日 报 ? 可 知 ” ,该 年 7 月 持续 高 
温 导 致 中 东部 受 旱 范围 和 程度 为 2015 一 2017 年 同 
期 之 最 ,影响 玉米 开花 授粉 ,导致 结实 率 下 降 , 同 时 
也 造成 马铃薯 结 暮 量 减 少 。2018 年 和 2019 年 河西 
地 区 的 干旱 程度 依然 严重 , 陇 中 、 陇 南 地 区 干旱 程 
度 较 2017 年 有 所 减轻 。 总 体 而 言 ,由 于 甘肃 省 复杂 
的 气候 类 型 广阔 的 面积 以 及 多 样 的 地 形 ,使 得 河 
西 地 区 常年 是 特 旱 区 和 重 旱 区 , 陇 中 地 区 是 重 旱 区 
和 中 晶 区 ,甘南 .陇南 和 陇 东 地 区 是 干旱 -无 旱 波 动 
变化 区 。 

2011 一 2019 年 ,甘肃 省 农业 干旱 程度 变化 在 农 
作物 生长 季 的 各 月 份 呈现 较 大 的 空间 异 质 性 (图 
6)。 以 CMDI 的 线性 倾向 率 作为 表征 指标 , 当 该 值 
为 负 时 ,表明 干旱 程度 增强 ;反之 意味 着 干旱 程度 
有 所 减轻 。 从 各 月 份 来 看 ,4 月 甘肃 省 大 部 分 区 域 
CMDI 的 线性 倾向 率 大 于 0, 说 明王 旱 程度 呈现 减轻 
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图 5 2011 一 2019 年 甘肃 省 农业 干旱 程度 
Fig. 3 Agricultural drought degree in Gansu Province during 2011—2019 


趋势 ,其 中 陇南 地 区 最 为 明显 ;5 月 除 陇 东 地 区 和 河 
西 地 区 西部 干旱 程度 略 有 减轻 外 ,河西 大 部 分 地 
区 、 陇 中 地 区 和 陇南 地 区 干旱 程度 有 所 加 剧 ;6 月 陇 
东 地 区 和 陇南 地 区 干旱 程度 明显 减轻 ,而 陇 中 地 区 
北部 和 河西 地 区 东部 的 石 羊 河流 域 下 游 旱 情 有 所 
增强 ;7 月 整个 河西 地 区 呈现 较 明显 的 干旱 程度 减 
轻 趋 势 , 陇 东 地 区 和 陇南 地 区 部 分 区 域 呈 现 干旱 程 
度 增 加 趋势 ;8 月 干旱 程度 明显 加 剧 的 区 域 为 陇南 
地 区 ,而 河西 地 区 的 干旱 程度 呈现 减轻 趋势 ;9 月 除 
陇南 和 陇 中 部 分 地 区 呈现 干旱 程度 减轻 趋势 外 ,其 
余地 区 均 呈 现 干旱 程度 增加 趋势 。 需 要 指出 的 是 ， 
4 一 9 月 甘南 地 区 干旱 程度 均 呈 减轻 趋势 。 


3.2.2 农业 干旱 时 间 变 化 规律 ”2011 一 2019 年 模拟 
CMDI 均 值 在 各 月 份 的 变化 较 大 (图 7) ,根据 值 的 大 
小 将 干旱 划分 为 特 旱 EE .中 旱 、 轻 旱 和 无 旱 5 个 
干旱 等 级 (图 8)。 如 图 7 所 示 ,CMDI 均 值 在 9 a 内 最 
大 值 为 1.10, 表 现 为 无 旱 , 则 表明 在 2014 年 7 月 无 任 
何 干 旱 胁迫 ;最 小 值 为 -2.54, 表 现 为 特 旱 , 则 表明 干 
旱 程度 在 2017 年 6 月 最 严重 。4 月 CMDI 均 值 介 
于 -1.25~0.33 之 间 , YE 2011,2018 年 表现 为 无 旱 ， 
2012 .2017 年 表现 为 中 旱 ,其 余年 份 均 表现 为 轻 旱 ; 
5 H CMDI 均 值 介 于 -1.82~ -0.06, 在 2017 年 表现 为 
重 旱 ,2012、2014、2019 年 表现 为 无 旱 ,其 余年 份 均 
表现 为 轻 旱 ;6 月 CMDI 均 值 介 于 -2.54~1.01, 在 
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图 6 2011 一 2019 年 4 一 9 月 甘肃 省 农业 干旱 变化 趋势 
Fig.6 Trends of agricultural drought in Gansu Province from April to September during 2011—2019 


2017 年 表现 为 特 旱 ,2019 年 表现 为 中 旱 2011, 
2013,2015 年 表现 为 轻 旱 ,其 余年 份 均 表现 为 无 旱 ; 
7 月 CMDI 均 值 介 于 -1.98~1.10, 在 2017 年 表现 为 重 
旱 ,2012 .2014 年 表现 为 无 旱 , 其 余年 份 均 表现 为 轻 
旱 ;8 月 CMDI 均 值 介 于 -1.22~1.05, 在 2017 年 表现 
为 中 旱 ,2011 .2014 .2016 年 表现 为 轻 旱 ,其 余年 份 
均 表 现 为 无 旱 ;9 月 均值 介 于 -1.32~1.03 ,在 2018 年 
表现 为 中 旱 ,2011、2017、2019 年 表现 为 轻 旱 ,其余 
年 份 均 表 现 为 无 旱 。 可 见 ,2012 年 的 干旱 程度 最 
轻 ,在 4 月 表现 为 中 旱 外 ,其 余 月 份 均 呈 现 无 旱 ; 
2017 年 的 干旱 程度 最 重 , 在 农作物 生长 季 (4 一 9 月 ) 
均 有 不 同 程度 的 干旱 发 生 。 


4 讨论 


通过 对 甘肃 省 农业 干旱 分 析 结 果 与 庞 素 非 
等 "研究 比较 发 现 ,甘肃 省 每 年 部 有 农业 干旱 发 
Æ ,河西 地 区 旱情 重 , 甘 南 地 区 旱情 轻 ,干旱 呈现 出 
明显 的 地 域 分 异 规律 , 即 “ 西 北 旱 .东南 润 ” ;本文 提 
取 甘 肃 省 农田 区 集成 多 因子 进行 监测 ,而 后 者 在 基 
于 单一 TVDI 指 数 的 基础 上 ,利用 相关 分 析 探 讨 海 


Jk 坡度 气温、 降水 .土地 利用 类 型 植被 类 型 与 农 
业 干 旱 的 密切 程度 ,本 研究 剔除 了 非 农 田 区 下 垫 面 
的 影响 ,避免 了 单一 遥感 干旱 因子 的 局 限 性 。 韩 兰 
英 等 利用 1960 一 2012 年 甘肃 省 农业 干旱 灾情 和 
气象 站 点 资料 ,构建 农业 干旱 灾害 风险 指数 ,以 作 
物 综合 损失 率 定量 分 析 农 业 干 旱 ,并 确定 气温 和 降 
水 是 关键 致 灾 因 子 。 与 之 相 比 ,本 文 模拟 的 CMDI 弥 
补 了 基于 气象 站 点 监测 不 连续 性 的 不 足 。 赖 力 等 " 
基于 SPEI 和 VCI 来 表征 气象 干旱 和 农业 干旱 ,发现 
农业 干旱 对 3 个 月 尺度 的 气象 干旱 反映 最 为 灵敏 。 
本 文 构 建 的 模型 中 耦合 了 VCI, TCI, PCI, TVDI, 
SMCI, DEM MIDI 等 因子 ,能 更 准确 地 反映 气象 干 
旱 是 导致 农业 干旱 的 外 在 驱动 力 ,揭示 农业 干旱 对 
气象 干旱 的 响应 。 同 时 ,针对 SPEI 与 模型 结果 的 相 
关 分 析 , 本文 仅 采 用 了 1 个 月 尺度 的 SPEI 值 , 若 采 
用 较 多 尺度 的 SPEI, 则 会 找 出 农业 干旱 对 气象 干旱 
反映 最 为 灵敏 的 时 间 尺 度 。 另 外 ,随机 森林 模型 的 
输入 参数 都 是 在 前 人 研究 的 基础 上 选取 的 干旱 因 
子 , 没 有 进行 设置 最 佳 模型 性 能 方案 ,因此 对 多 因 
子 进 行 合 理 的 筛选 验证 ,利用 更 多 成 熟 的 遥感 致 旱 
子 输入 随机 森林 模型 估算 农业 干旱 指数 是 今后 
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图 7 2011 一 2019 年 4 一 9 月 模拟 CMDI 均 值 


Fig.7 Averaged simulated CMDI from April to September during 2011—2019 


图 8 不 同年 月 干旱 程度 


Fig.8 Drought degree in different years or months 


需要 考虑 的 方向 。 
甘肃 省 河西 地 区 是 农业 干旱 发 生 的 特 旱 区 ， 
FRANE .土地 盐 碱 化 ,水 土 流失 严重 等 问题 使 得 


绿洲 成 为 维系 人 类 生产 生活 的 环境 基础 。 若 制定 
合理 的 生态 保护 政策 .限制 绿洲 资源 的 过 度 开 发 、 
广泛 发 展 戈壁 农业 和 设施 农业 等 一 系列 节 水 措施 ， 
则 可 缓解 农业 干旱 ;从 长 期 来 看 ,还 可 增加 该 地 区 
耐 旱 作物 的 种 植 面积 和 实现 跨 流 域 调 水 工程 这 一 
目标 减缓 农业 和 干旱 。 陇 中 地 区 是 农业 干旱 发 生 的 
重 旱 区 (包括 庆 阳 市 .平凉 市 .定西 市 ) ,由 于 降雨 量 
^H ER .地 下 水 资源 不 多 、 基 本 没有 外 来 水 源 和 大 型 
河流 等 因素 限制 ,使 得 该 区 域 的 水 资源 严重 匮乏 。 
若 能 高 效 利 用 黄河 上 游 流域 水 源 、 加 强 流 域 生态 保 
JP RERA KEREDE .广泛 发 展 温 室 农业 和 灌 
溉 节 水 农业 , 则 可 以 有 效 减缓 该 地 区 的 农业 干旱 。 


5 结论 

本 文 利 用 多 源 遥 感 数 据 提取 农业 致 旱 因 子 , 结 
合 随机 森林 模型 构建 CMDI, 并 以 实测 SPEI 和 CMDI 
来 评价 模型 结果 ,最 后 分 析 甘 肃 省 农业 干旱 的 空间 
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分 异 特征 和 变化 规律 ,得 到 以 下 主要 结 

(1) 模型 结果 与 实测 A ERTE 
0.935 ,最 低 为 0.634, 精 度 较 高 ;与 SPEI 相 关 性 大 多 
在 0.50 以 上 ,空间 上 可 靠 性 强 。 这 表明 利用 随机 和 森 
林 模 型 集成 多 因子 监测 农业 干旱 是 有 效 的 ,可 以 很 
好 地 反映 农业 干旱 的 发 生发 展 过 程 。 

(2) 2011 一 2019 年 ,甘肃 省 农业 干旱 程度 最 为 
严重 的 年 份 是 2017 年 ,最 轻 的 年 份 是 2012 年 ;甘肃 
省 农业 干旱 呈现 较 大 的 波动 趋势 ,年 际 变化 较为 复 
2 ,但 呈现 出 明显 的 地 域 分 异 规 律 , 即 由 甘肃 省 东 
南部 向 西北 部 干旱 程度 加 重 。 其 中 ,河西 地 区 多 为 
特 旱 区 和 重 旱 区 , 陇 中 地 区 为 重 (中 ) 旱 区 ,陇南 、 陇 
AR 上 冉 南 地 区 为 干旱 -无 早 波动 变化 区 。 

(3) 2011—2019 年 ,甘南 地 区 在 农作物 生长 季 
干旱 程度 均 呈 减轻 趋势 , 除 大 部 分 地 区 在 4 月 和 6 
月 ,以 及 陇 东 地 区 和 陇南 地 区 分 别 在 5 月 和 9 月 农 
业 旱 情 有 所 减轻 趋势 外 ,其 余地 区 在 农作物 生长 季 
的 旱情 呈 增 加 趋势 。 


参考 文献 (References) 


[1] KIER, 朱 秀 芳 , 刘 婷 婷 , 等 . 基于 多 维 Copula 的 中 国 干旱 特征 
及 危险 性 分 析 町 . 干旱 区 地 理 , 2022, 45(2): 333-345. [Zhang 
Shizhe, Zhu Xiufang, Liu Tingting, et al. Drought characteristics 


and risk hazard in China based on multidimensional Copula model 

[J]. Arid Land Geography, 2022, 45(2): 333-345. ] 

nun 青海 省 干旱 时 空 特征 及 与 大 气 环流 响应 
资源 与 环境 , 2021, 35(12): 60-65. [Fan Lei, Lü 


p) Jis HAM 
关系 四. 干旱 


Aifeng, Zhang Wenxiang. Temporal- spatial variation characteris- 


tics of drought and its relationship with atmospheric circulation in 
Qinghai Province[J]. Journal of Arid Land Resources and Environ- 
ment, 2021, 35(12): 60-65. | 

[3] IPCC. Climate Change 2021: The physical science basis: Contribu- 
tion of working group I to the sixth assessment report of the inter- 
governmental panel on climate change[M]. Cambridge: Cambridge 
University Press, 2021. 

[4]  Thilagaraj P, Masilamani P, Venkatesh R, et al. Google earth en- 
gine based agricultural drought monitoring in Kodavanar Water- 
shed, part of Amaravathi Basin, Tamil Nadu, India[J]. Internation- 
al Archives of the Photogrammetry, Remote Sensing and Spatial 
Information Sciences, 2021, XLIII- B5-2021: 43-49. 

[5] FNA, 陈云 浩 , 孙 洪 果 . 典型 农业 干旱 遥感 监测 指数 的 比较 及 分 
35 Ms ZR UU]. 农业 工程 学 报 , 2012, 28(14): 147-154. [Sun Hao, 


Chen Yunhao, Sun Hongquan. Comparisons and classification sys- 


tem of typical remote sensing indexes for agricultural drought[J]. 


Transactions of the Chinese Society of Agricultural Engineering, 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


2012, 28(14): 147-154. ] 
吴 泽 棉 , 印 建 秀 , 刘 苏 峡 , 等 . 基于 土壤 水 分 的 农业 干旱 监测 研 
SEE]. 地 理科 学 进展 , 2020, 39(10): 1758-1769. [Wu Zemi- 


an, Qiu Jianxiu, Liu Suxia, et al. Advances in agricultural drought 


monitoring based on soil moisture[lJ]. Progress in Geography, 
2020, 39(10): 1758-1769. ] 

Dai A. Increasing drought under global warming in observations 
and models[J]. Nature Climate Change, 2013, 3(1): 52-58. 

刘 完 锋 , REI, TIR, 等 . 农业 干旱 监测 研究 进展 与 展望 吊 . 
地 理学 报 , 2015, 70(11): 1835-1848. [Liu Xianfeng, Zhu Xiufang, 


Pan Yaozhong, et al. Agricultural drought monitor: Progress, chal- 
lenges and prospect[J]. Acta Geographica Sinica, 2015, 70(11): 
1835-1848. ] 

Lee S J, Kim N R, Lee Y W. Development of integrated crop 
drought index by combining rainfall, land surface temperature, 
evapotranspiration, soil moisture, and vegetation index for agricul- 
tural drought monitoring[J]. Remote Sensing, 2021, 13(9): 1778, 
doi: 10.3390/1513091778. 

WIE K, 李 净 , 王 丹 , 等 . 基于 MODIS 和 TRMM 数据 的 黄土 高 原 
农业 干旱 监测 四. 干旱 区 < 地 更 2019, 42(1): 173-178. [Hu Peng- 
fei, Li Jing, Wang Dan, et al. Monitoring agricultural drought in 
the Loess Plateau using MODIS and TRMM data[J]. Arid Land Ge- 
ography, 2019, 42(1): 173-178. ] 

温 庆 志 , PMG, 张强 , 等 . 基于 多 源 遥 感 数据 的 农业 干旱 监测 模 
型 构建 及 应 用 员 . 生态 学 报 , 2019, 39(20): 7757-7770. [Wen 
Qingzhi, Sun Peng, Zhang Qiang, et al. An integrated agricultural 


drought monitoring model based on multi-source remote sensing 
data: Model development and application[J]. Acta Ecologica Sini- 
ca, 2019, 39(20): 7757-7770. ] 

刘 高 鸣 , 谢 传 节 , 何 天 乐 . 基于 多 源 数据 的 农业 干旱 监测 模型 构 
EJ]. 地 球 信 息 科 学 学 报 , 2019, 21(11): 1811-1822. [Liu Gaom- 


ing, Xie Chuanjie, He Tianle. Agricultural drought monitoring 


model constructing based on multi-source data[J]. Journal of Geo- 
information Science, 2019, 21(11): 1811-1822. ] 

Brown J F, Wardlow B D, Tadesse T, et al. The vegetation drought 
response index (VegDRI): A new integrated approach for monitor- 
ing drought stress in vegetation[JJ. GlScience & Remote Sensing, 
2008, 45(1): 16-46. 

Wu J, Zhou L, Liu M, et al. Establishing and assessing the integrat- 
ed surface drought index (ISDI) for agricultural drought monitoring 
in mid-eastern China[J]. International Journal of Applied Earth 
Observation and Geoinformation, 2013, 23: 397-410. 

EBER, 李 晓 松 , VERS, 等 . 光学 与 雷达 遥感 协同 的 大 尺度 草地 
灌 从 化 监测 研究 四 . 干旱 区 资源 与 环境 , 2021, 35(2): 130-135. 


[Shao Jing, Li Xiaosong, Yang Junting, et al. Study on large scale 


grassland shrub monitoring based on optical and radar remote 
sensing[J]. Journal of Arid Land Resources and Environment, 
2021, 35(2): 130-135. ] 

刘斌 , WE, RFE. 基于 随机 森林 模型 的 土壤 


iml 


E 金 属 源 解 析 


202302.00276v1 


chinaXiv 


[17] 


[18] 


[19] 


[20] 


21] 


[22] 


3] 李 


[24] 


[25] 


8] K BT: METER 


— UE PRAA. 干旱 区 资源 与 环境 , 2019, 33(1): 106- 
111. [Liu Bin, Guo Xing, Zhu Yu’ en. Analysis of soil heavy metal 


sources in Jinzhong Basin based on random forest model[J]. Journal 
of Arid Land Resources and Environment, 2019, 33(1): 106-111. ] 

高 亮 , FE. 李强 子 , 等 . 融合 土地 覆盖 和 土壤 水 分 产品 的 近 地 
表 空 气温 度 空 间 化 方法 中 . 地 球 信息 科学 学 报 , 2020, 22(10): 
2023-2037. [Gao Liang, Du Xin, Li Qiangzi, et al. A near-surface 


air temperature spatialization method integrating landuse and soil 
moisture products[J]. Journal of Geo-information Science, 2020, 22 
(10): 2023-2037. ] 

沈 润 平 , 郭 佳 , IKAR. 基于 随机 森林 的 遥感 干旱 监测 模型 的 构 
Æ. 地 球 信 息 科 学 学 报 , 2016, 19(1): 125-133. [Shen Run- 


ping, Guo Jia, Zhang Jingxian. Construction of a drought monitor- 


ing model using the random forest based remote sensing[J]. Jour- 
nal of Geo-information Science, 2016, 19(1): 125-133. ] 

董 婷 , 任 东 , 邵 攀 , 等 . 基于 多 源 遥 感 数据 和 随机 森林 的 综合 旱 
TETE bs F3 ELT]. 农业 机 械 学 报 , 2019, 50(8): 200-212. [Dong 


Ting, Ren Dong, Shao Pan, et al. Construction of integrated 


drought condition index based on multi-sensor remote sensing and 

random forest[J]. Transactions of the Chinese Society for Agricul- 

tural Machinery, 2019, 50(8): 200-212. | 

VERE, 魏 伟 , 郭 泽 呈 , 等 . 基于 TVDI 的 甘肃 省 农业 旱情 特征 及 
其 影响 因素 中. adore. 2019, 38(6): 1849-1860. [Pang Suf- 


ei, Wei Wei, Guo Zecheng, et al. Agricultural drought characteris- 


tics and its influencing factors in Gansu Province based on TVDI 
[J]. Chinese Journal of Ecology, 2019, 38(6): 1849-1860. ] 

胡蝶 , vb, 王丽娟 , 等 . 欧 空 局 主 被 动 微波 土壤 水 分 融合 产品 
在 甘肃 省 干旱 监测 中 的 应 用 町 . 干旱 气象 , 2019, 37(4): 517- 
528. [Hu Die, Sha Sha, Wang Lijuan, et al. Drought monitoring ap- 


plications of the European space agency climate change initiative 
soil moisture combined product in Gansu Province[J]. Journal of 
Arid Meteorology, 2019, 37(4): 517—528. | 
HA, 张 勃 , 马 尚 谦 , 等 . 甘肃 省 河东 地 区 春玉 米 干旱 演变 特征 
及 灾 损 风险 区 划 町 . 干旱 地 区 农业 研究 , 2019, 37(6): 238-247. 


[Wei Kun, Zhang Bo, Ma Shangqian, et al. Drought evolution char- 


acteristics and risk zoning of spring maize in Hedong of Gansu 
Province[J]. Agricultural Research in the Arid Areas, 2019, 37(6): 
238-247. ] 

:栋梁 , 刘 德 祥 . 甘肃 气候 [M]. 北京 : 气象 出 版 社 , 2000: 2-105. 


[Li Dongliang, Liu Dexiang. Gansu climate[M]. Beijing: Meteoro- 
logical Publishing House, 2000: 2-105. | 

郭 斌 , 邢 文 雪 , 李 云 , 等 . 基于 TRMM 3B43 产品 的 1998—2016 
年 黄 淮 海平 原 降水 时 空 特征 研究 四. 地 理科 学 , 2021, 41(7): 
1285-1294. [Guo Bin, Xing Wenxue, Li Yun, et al. Temporal and 


spatial characteristics of precipitation based on TRMM 3B43 prod- 
uct in the Huang- Huai- Hai Plain during 1998—20196][J]. Scientia 
Geographica Sinica, 2021, 41(7): 1285-1294. ] 

HIRE, 李 艳 忠 , 余 文 君 , 等 . 遥感 降水 产品 在 中 国 不 同 气候 区 
的 适用 性 研究 四. 地 球 信 息 科 学 学 报 , 2021, 23(7): 1296-1311. 
[Peng Zhenhua, Li Yanzhong, Yu Wenjun, et al. Research on the 


[5 
CA 


数据 的 甘肃 省 农业 干 


[26] 


[27] 


[28] 


[29] 


[30] 


[31] 


[32] 


[33] 


[34] 


[35] 


ChinaXiv 合 作 期 刊 


旱 研究 


21 


出 


applicability of remote sensing precipitation products in different 
climatic regions of China[J]J. Journal of Geo- information Science, 
2021, 23(7): 1296-1311. ] 
费 龙 , 邓 国 荣 ， jc 
时 空 特征 四. 自然 资 


Long, Deng Guorong, Zhang Hongyan, et al. Spatiotemporal pat- 


E. AE TREK Z TEC SUL EE EEUU 5:83 
^K JR, 2020, 35(12): 3051-3063. [Fei 


terns of precipitation and drought and flood using Z-index in Dem- 
ocratic People’s Republic of Korea[J]. Journal of Natural Resourc- 
es, 2020, 35(12): 3051-3063. | 

ARCET, 李 丽 娟 , 李 九 一 . 基于 TRMM 降 尺度 和 MODIS 数 据 的 
综合 干旱 监测 模型 构建 中 . 自然 资源 学 报 , 2020, 35(10): 2553- 


2568. [Yu Haozhe, Li Lijuan, Li Jiuyi. Establishment of compre- 


hensive drought monitoring model based on downscaling TRMM 
and MODIS data[J]. Journal of Natural Resources, 2020, 35(10): 
2553-2568. | 

李 净 , 张 晓 . TRMM [ez ts BI) 7 ENEIT ADAI]. 地 理科 
学 , 2015, 35(9): 1164-1169. [Li Jing, Zhang Xiao. Downscaling 
method of TRMM satellite precipitation data[J]. Scientia Geograph- 
ica Sinica, 2015, 35(9): 1164-1169. ] 

Yu G, Di L. Downscaling of global soil moisture using auxiliary da- 
ta[ C//Yang W. International Geoscience & Remote Sensing Sym- 
posium. Institute of Electrical and Electronics Engineers, 2008: 
230-233. 

杜 灵通 , MKA, Efi, 等. AET Ae ot due t 
Ji PEEL. 农业 工程 学 报 , 2014, 30(9): 126-132. [Du Ling- 


tong, Tian Qingjiu, Wang Lei, et al. A synthesized drought moni- 


E Ws 监测 


toring model based on multi-source remote sensing data[J]. Trans- 

actions of the Chinese Society of Agricultural Engineering, 2014, 

30(9): 126-132. ] 

IKER, 魏 伟 , DES E, 等 . 基于 FY-3C 和 TRMM 数据 的 西北 干旱 
X F Uu n 5; PE [I]. 生态 学 杂志 , 2020, 39(2): 690-702. 

Zhang Jing, Wei Wei, Pang Sufei, et al. Monitoring and assess- 


ment of drought in arid area in northwest China based on FY-3C 
and TRMM Data[J]. Chinese Journal of Ecology, 2020, 39(2): 690- 
702. ] 
谢 五 三 , 田 红 , E JE. 基于 CI 指数 的 淮河 流域 干旱 时 空 特征 研 
5E [J]. 气象 , 2013, 39(9): 1171-1175. [Xie Wusan, Tian Hong, 


Wang Sheng. Study on temporal-spatial characteristics of drought 
in Huaihe River Basin based on CI index[J]，Meteorological 
Monthly, 2013, 39(9): 1171-1175. ] 

GBT 20481. 2006 气象 干旱 等 级 [S]. 2006: 40-41. [GBT 20481. 
Meteorological drought grade in 2006[S]. 2006: 40-41. ] 

张 调 风 , IKE, 王 有 恒 , 等 . 基于 综合 气象 干旱 指数 的 石 羊 河流 
域 近 50 年 气象 干旱 特征 分 析 叫 . 生态 学 报 , 2013, 33(3): 975- 
984. [Zhang Tiaofeng, Zhang Bo, Wang Youheng, et al. Drought 


characteristics in the Shiyang River Basin during the recent 50 
years based on a composite index[J]. Acta Ecologica Sinica, 2013, 
33(3): 975-984. | 

Vicente-Serrano S M, Beguería S, López-Moreno J I. A multiscalar 


drought index sensitive to global warming: The standardized pre- 


202302.00276v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


22 tib E! 4636 
cipitationevapotranspirationindex-SPEI[J]. Journal of Climate, 2010, 2018- 01- 04(008). [Hong Wenquan. Gansu released the top 10 
23(7): 1696-1718. weather and climate events of 2017[N]. Gansu Daily, 2018-01-04 

[36] 27s, 温 兆 飞 , 陈 吉 龙 , 等 . 基于 随机 森林 回归 模型 的 悬浮 泥 (008). ] 
沙 浓度 遥感 估算 册 . 遥感 学 报 , 2019, 23(4): 756-772. [Fang Xin- [39] HESE, 张强 , 赵 红 岩 , 等 . 甘肃 省 农业 干旱 灾害 损失 特征 及 其 
rui, Wen Zhaofei, Chen Jilong, et al. Remote sensing estimation of 对 气候 变 暧 的 响应 四 . 中 国 沙漠 , 2016, 36(3): 767-776. [Han 
suspended sediment concentration based on random forest regres- Lanying, Zhang Qiang, Zhao Hongyan, et al. The characteristics 
sion modellJ]. National Remote Sensing Bulletin, 2019, 23(4): agricultural drought disaster loss and response to climate warming 
756-772. | in Gansu, China[J]. Journal of Desert Research, 2016, 36(3): 767- 
[37] ftm, 杨 联 安 , EWE, 等 . 基于 随机 森林 的 苹果 区 土壤 有 机 质 空 716.] 
[RESCUE]. 干旱 区 资源 与 环境 , 2018, 32(8): 141-146. [Ren Li, [40] HJ, EE, 冯 凯 . 甘肃 省 农业 干旱 对 多 尺度 气象 干旱 的 响应 
Yang Lian’ an, Wang Hui, et al. Spatial prediction of soil organic [J]. 节 水 灌溉 , 2020(2): 102-108. [Lai Li, Su Xiaoling, Feng Kai. 
matter in apple region based on random forest[JJ. Journal of Arid Response of agricultural drought to multi-scale meteorological 
Land Resources and Environment, 2018, 32(8): 141-146. | drought in Gansu Province[J]. Water Saving Irrigation, 2020(2): 


| 


[38] WR. 甘肃 发 布 2017 年 十 大 天 气 气候 事件 [N]. 甘肃 日 报 ， 102-108. | 


Agricultural drought research based on multi-source remote sensing data in 
Gansu Province 


ZHANG Juan, YAO Xiaojun LI Jing WANG Xiaoyan 


(College of Geography and Environmental Science, Northwest Normal University, Lanzhou 730070, Gansu, China) 


Abstract: Drought is the main environmental stressor to crop growth and development and a pivotal natural 
constraint to high agricultural yield and harvest. Agricultural drought monitoring is typically based on 
meteorological observations, which, to some extent, cannot reflect the agricultural drought conditions on a 
regional scale. Taking the Gansu Province of China as the study area, the compositive meteorological drought 
index (CMDI) was constructed using a random forest regression model and the spatiotemporal pattern and change 
rule of drought during the crop-growing season (April to September) from 2011 to 2019 in the Gansu Province 
were analyzed based on MODIS, TRMM, ESA CCI, and other remote sensing data products as well as 
meteorological station data. The results are as follows: (1) the coefficient of determination (R^) of the CMDI and 
the measured value were all greater than 0.634 in each month; the CMDI has a high spatial correlation with the 
standardized precipitation evapotranspiration index, showing that the CMDI can reflect the occurrence and 
development of agricultural drought. (2) The agricultural drought in the Gansu Province exhibited an obvious 
pattern of regional differentiation, with the drought degree gradually increasing from southeast to northwest. 
Most parts of the Hexi region were in special and severe drought, the Longzhong region was in severe (medium) 
drought, and the Longdong, Longnan, and Gannan regions were in fluctuating drought. (3) From 2011 to 2019, 
the agricultural drought in the Gansu Province exhibited a significant fluctuating trend on both annual and 
monthly scales, and the drought degree was most moderate in 2012 and most severe in 2017. Agricultural drought 
has been reduced in the Gannan region in April-September, most areas in April and June, the Longdong region in 
May, and the Longnan region in September. However, the remaining regions experienced an increasing trend of 
drought during the crop-growing season from 2011 to 2019. 
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